Magnon dispersion and thermodynamics in CsNiF 3 
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We present an accurate transfer matrix renormalization group calculation of the thermodynamics 
in a quantum spin-1 planar ferromagnetic chain. We also calculate the field dependence of the 
magnon gap and confirm the accuracy of the magnon dispersion derived earlier through an 1/n 
expansion. We are thus able to examine the validity of a number of previous calculations and further 
analyze a wide range of experiments on CsNiF3 concerning the magnon dispersion, magnetization, 
susceptibility, and specific heat. Although it is not possible to account for all data with a single 
set of parameters, the overall qualitative agreement is good and the remaining discrepancies may 
reflect departure from ideal quasi-one-dimensional model behavior. Finally, we present some indirect 
evidence to the effect that the popular interpretation of the excess specific heat in terms of sine- 
Gordon solitons may not be appropriate. 
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I. INTRODUCTION 

The magnetic compound CSN1F3 undergoes three- 
dimensional (3D) ordering at very low temperatures, 
T < Tn = 2.7 K, but exhibits essentially one-dimensional 
(ID) behavioc. for T > Tn- A number of experimental 
investigationsa suggest that an appropriate ID model is 
described by the spin s = 1 Hamiltonian 



W 



[-JS n -S n+1 + A(S^) 2 - 9 /x B H-S n ] , (1.1) 



which contains a ferromagnetic (J > 0) isotropic ex- 
change interaction and an easy-plane (A > 0) single-ion 
anisotropy, in addition to the usual Zeeman term pro- 
duced by an applied field H. 

The derivation of accurate theoretical predictions 



based on Hamiltonian (1.1) turned out to be more diffi- 
cult than anticipated thanks to the strong quantum fluc- 
tuations that occur in this quasi-lD system. In particu- 
lar, the leading-order magnon dispersion derived within 
the usual 1/s expansion is too crude an approximation 
for s = 1. As a result, inelastic neutron scattering ex- 
periments were analyzeda mostly on the basis of an ak 
ternative dispersion derived by Lindgard and Kowalskad 
using a self-consistent approach that is designed to prop- 
erly account for single-ion anisotropy. Similarly, a large 
body of experimental data became available for thermo- 
dynamic quantities such as magnetization, susceptibility, 
and specific heat, but a corresponding theoretical calcula- 
tion proceeded slowly. To the best of our knowledge, the 
most accurate calculation of thermodynamics was pro- 
vided by Delica et al.cl based on a quantum transfer ma- 
trix, while comparable success was claimed more recently 
by Cuccoli et allil through a sophisticated semiclassical 
approach. The above two papers also contain an exten- 
sive list of references to earlier work. 



It is the aim of the present paper to derive theoretical 
predictions that are accurate to within line thickness and 
thus provide a safe basis for the discussion of various 
issues that have been raised during the long history of 
this subject. 

In Sec. Q, experimental data on the magnon disper- 
sion are analyzed in terms of an unconventional 1/n 
expansionQ which is shown to contain the Lindgard- 
Kowalska dispersion as a special case. The accuracy 
of the leading 1/n approximation is confirmed by an 
independent calculation of the field dependence of the 
magnon gap using a density matrix renormalization 
group (DMRG) methodjll while a discussion of anhar- 
monic corrections within the conventional 1/s expansion 
is also included for comparison. Thermodynamic quanti- 



ties are calculated in Sec. [II by a powerful transfer ma- 
trix renormalization group (TMRG) algorithmOil3 which 
addresses directly the infinite-chain limit. We are thus in 
a position to appreciate the relative accuracy of earlier 
calculations, analyze all available data and anticipate re- 
sults of possible future experiments, as well as challenge 
popular interpretations in terms of sine-Gordon solitons. 
A brief summary of the main conclusions is given in Sec. 
[V. 



II. THE MAGNON DISPERSION 

The standard spinwave theory is a method for calcu- 
lating quantum corr ections around the classical minimum 
of Hamiltonian (1.1) by a systematic 1/s expansion. The 
1 /n expansion developed in Ref . ^| is of a similar nature, 
except that the corresponding "classical" minimum is a 
variational Hartree-like ground state that is more sensi- 
tive to the nature of single-ion anisotropy and thus pro- 
vides a more sensible starting point. Hence one obtains 
an accurate magnon dispersion even if the 1/n series is 



1 



restricted to the harmonic approximation. 

For a field applied in a direction perpendicular to the 
c-axis, e.g., H = (H, 0,0), the magnon energy at crystal 
momentum q is given by 



2J 



{(l+e)(J-cos g ) 



;) cos q I 
(2.1) 

Here and in the rest of the paper we employ rationalized 
parameters for anisotropy and field, 



1/2 



= A/J, h = gj_n B H/ J, 



(2.2) 



while energy and temperature may be measured in units 
of the exchange constant J. The notation employed for 
the gyromagnetic ratio g± implies that the corresponding 
ratio <7|i for a field parallel to the c-axis may be di ffere nt. 



Finally, the dimensionless parameter e in Eq. ( 2.1 



determined in terms of a and h by the algebraic equation 



q(l-e 2 ) 1/2 

2/l + 4(l-£ 2 ) 1 /2' 



(2.3) 



One should add that derivation of systematic 1 jn correc- 
tions to the harmonic approximation (2.1) is possible^ 
but unnecessary in the parameter range of current inter- 
est: a,h < 0.5. 

At zero field, the root of Eq. (2.3) is e = a/4 which 
is inserted in Eq. (2.1) to provide a completely explicit 
expr essio n for the magnon dispersion. For nonzero field, 
Eq. (2.3) may be solved by simple iteration starting with 
e = 0. In fact, the result of a single iteration, 

(2.4) 



2h + 4 ' 

i s pr actically indistinguishable from the exact root of Eq. 

(2.3) for parameters such that a,h < 0.5. The last re- 
mark becomes especially important if one notes that the 
disp ersion obt aine d by inserting the approximate root 

(2.4) in Eq. (2.1) is precisely the magnon dispersion 
derived earlier by Lindgard and Kowalskajj applied for 
s = 1, which was in turn employed for the analysis of 



experimental data from inelastic neutron scattering! 

The latter analysis provided what is often referred to 
as the standard set of parameters for CSN1F3: 



J = 23.6K, A = 9K, 



= 2.4. 



(2.5) 



The corresponding theoretical predictions of the mag-aon 
dispersion (2.1) are compared to experimental dataa in 
the upper panel of Fig. EL The agreement is obviously 
very good for field H = 41 kG, while a slight but system- 
atic deviation is observed for H = 0. This conclusion is 
somewhat surprising in view of the claim in Ref. |^ that 
nearly perfect agreement is obtained for both field values, 
even though the Lindgard-Kowalska dispersion empl oyed 
in the above reference is practically identical to Eq. (2.1) 




FIG. 1. The magnon energy E = u) n a s a function of 
crystal momentum q calculated from Eq. (2.1) for two values 
of the applied field, H — and H — 41 kG, and two different 
sets of parameters. The insets illustrate the corresponding 
field dep endence of the q = magnon gap G calculated from 
Eq. (2.7). Solid circles represent experimental data from Ref. 
I taken at T = 4.2 K. 



for the set of parameters (2.5). The systematic nature of 
this discrepancy makes it-junlikely that the data commu- 
nicated to us by Steinerliil differs from the data actually 
used in the analysis of Ref. || A more likely explanation 
is that the Lindgard-Kowalska dispersion was further ap- 
proximated by the authors of Ref. |[ as is evident in the 
expression for the magnon gap given in their Eq. (5). 

Although the observed discrepancy appears to be mi- 
nor, it nonetheless leads to a substantial redefinition of 
parameters. Thus we have redetermined the exchange 
constant J and anisotropy A by a least-square fit of the 
zero-field data to dispersion (2.1), while the gyromag- 
netic ratio was subsequently obtained by a one-parameter 
least-square fit of the H = 41 kG data. The resulting new 
set of parameters 



J = 23.6 K, A = 11 K, 



2.18 



(2.6) 



restores agreement with experiment for both field values, 
as is shown i n th e low er p anel of Fig. A notable fea- 
ture of Eqs. (2.5) and (2.6) is that the exchange constant 
has remained unchanged. Indeed, throughout our analy- 
sis, we found no evidence for departure of the exchange 
constant from the value J = 23.6 K which will thus be 
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adopted in the following without further questioning. 

In contrast, the observed significant fluctuations in the 
anisotropy constant A and gyromagnetic ratio g± sim- 
ply reflect the fact that the magnon dispersion is not 
especially sensitive to thos e pa ram eter s. Therefore, their 
values given in either Eq. (2.5) or (2.6) cannot be consid- 
ered as established without further corroboration. Now, 
the reduced value of the gyromagnetic ratio given in Eq. 
( |2.6D is consistent with g± = 2.1 ± 0.05 obtained inde- 
pendently by-, measuring the saturation magnetization at 
strong fieldsQ and is also supported by the analysis of 
the zero-field susceptibility in Sec. III. But a proper 



choice of the anisotropy constant A will be a matter 
of debate throughout this paper. In this respect, one 
should keep in mind that the neutron data displayed 
in Fig. [j] were taken at helium temperature, T = 4.2 
K, which is relatively high but not too distant from the 
3D-ordering transition temperature T/v — 2.7 K. Hence, 
finite-temperature effects as well as deviations from ideal 
ID behavior may already be present. 

An important special case of the magnon dispersion 
(|2.lD is the zero- momentum gap G = oj q= o, or 



g±[i B H + A 



1-e 



1/2' 



1/2 



(2.7) 



where we have made use of the algebraic equation ( |2.2| ) 
to simplify the expression^ A comparison of the predic- 
tions of Eq. (2/7) rjfi^h the measured field dependence 
of the magnon gapallil is shown in the insets of Fig. |l| 
for both sets of parameters. Although the overall agree- 
ment is reasonable, systematic deviations are present at 
relatively low field values in both cases. An attempt to 
redetermine the pa ram eters by a least-square fit of the 
q = data to Eq. (2.7) yields values for A and g± that 
would significantly compromise the agreement obtained 
at nonzero crystal momentum q. 

Implicit in the preceding discuss ion is the presump- 
tion that the magnon dispersion (2.1) and its special 
case (2.7) are sufficiently accurate and there is no need 
to proceed with the calculation of anharmonic 1/n cor- 
rections. We now test this assumption by a completely 
independent calculation of the field dependence of the 
magnon gap based on a density matrix reuarmalization 
group (DMRG) algorithmJZl An early effortll3 to apply a 
renormalization-group technique was restricted to short 
chains (16 sites) and thus provided reasonable but not 
especially accurate estimates of the magnon gap. The 
DMRG algorithm allowed us to calculate the gap on 
long chains up to 400 sites. We have also tested the 
stability of our results through Shanks or Richardson 
extrapolation^ and believe to have calculated the gap 
to an accuracy greater than the three figures actually 
displayed in the third column of Table |. 

It is then important that the corresponding results ob- 
tained through Eq. (2.7), listed in the second column 
of Table |, are in agreement with the DMRG calculation. 



TABLE I. Magnon gap in units of J, for a typical 
anisotropy a — A/ J = 0.38, and a field h — g±/j,BH/J ap- 
plied in a direction perpendicular to the c-axis. 



n 


1/n 


Magnon gap G 
DMRG 


1/s 


0.000 


0.000 


0.000 


0.000 


0.025 


0.105 


0.106 


0.109 


0.050 


0.152 


0.155 


0.160 


0.075 


0.192 


0.195 


0.201 


0.100 


0.227 


0.230 


0.238 


0.150 


0.290 


0.295 


0.304 


0.200 


0.350 


0.354 


0.365 


0.250 


0.406 


0.411 


0.422 


0.300 


0.461 


0.466 


0.478 


0.400 


0.568 


0.573 


0.586 


0.500 


0.673 


0.677 


0.691 



Since the relative accuracy is expected to further improve 
at nonzero crystal mom ent um q, one must conclude that 
the magnon dispersion (2.1) is sufficiently accurate for all 
practical purposes. Therefore, any disagreement between 
theory and experiment should be attributed to other rea- 
sons. In particular, one should note in Table || that the 
1/n results slightly underestimate the DMRG data and 
hence the latter cannot be used to eliminate the remain- 
ing small disagreement with the experimental data shown 
in the insets of Fig. [j]. 

Next we comment on the relative validity of the stan- 
dard semiclassical theory based on a 1/s expansion. The 
corresponding harmonic approximation of the magnon 
dispersion is clearly inaccurate, as is apparent in the 
estimate of anisotropy A — 4.5 K encountered in the 
early literature. El However, the semiclassical prediction 
can be significantly improved by including the first (an- 
harmonic) 1/s correction. At zero field, a completely 
analytical calculation is possible and may be found in 
Ref. For nonzero field, the anharmonic correction is 
expressed in terms of complicated integrals that cannot 
be computed analytically. Therefore, for simplicity, the 
main point is made here by considering only the q = 
magnon gap which can be written as 



G = Go [1 + S/s + 0(1/ s 2 )] 
G = sJ [h(h+ 2a)} 1 12 , 



If, 1 -cosq + h/2 + a/A 

I=~ dq -p-, 

n Jo [(1 -cosq + h/2)(l -cosq + h/2 + a)} ' 

where the rationalized field is now defined as h = 
OjU bH/sJ which differs from the definition given in Eq. 
( |2.2| ) by a factor that becomes unimportant for s = 1 . Go 
is the (harmonic) classical approximation and 6 provides 
the first anharmonic correction which amounts to about 
15 ~ 20% of the total ans wer. Numerical values for the 
gap calculated from Eq. (2.8), applied for s = 1, are 
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listed in the fourth column of Table |. These values over- 
estimate the DMRG data by a wider margin than the 
harmonic 1/n approximation underestimates the same 
data. T here fore, we again conclude tha t the magnon dis- 
persion (2.1) and the magnon gap ( [2.7| ) provide the most 
accurate description. 

Finally, we mention that an 1 /n expansion is also pos- 
sible in the case of a field parallel to the c-axis, along the 
lines outlined in the Appendix of Ref. [| Such a possi- 
bility will not be pursued further in the present paper, 
except for a minor application in Sec. |III B| , mainly be- 
cause we do not know of an experimental measurement 
of the magnon dispersion for this field orientation. 



III. THERMODYNAMICS 

The most straightforward method for calculating the 
partition function is a complete numerical diagonaliza- 
tion of the Hamiltonian on finite chains. The size of the 
resulting matrices is 3 N x 3 N and grows exponentially 
with the total number of sites N. Therefore, a calcu- 
lation is possible only on short chains while a reliable 
extrapolation to larger values of N is difficult. 

More powerful numerical methods proceed with the 
construction of a quantum transfer matrix (QTM) ob- 
tained by an M-step Trotter decomposition. An explicit 
calculation waSrijiitially performed via Quantum Monte 
Carlo samplingeJ and was also limited to short chains 
(N = 16) and a relatively small number of Trotter steps 
(M = 12). This procedure led to reasonable results for 
the magnetization and susceptibility, but the calculation 
of the specific heat was plagued by large statistical errors. 

A more systematic QTM calculation was later 
accomplisheda on long chains (N ~ 150) by limiting 
the number of Trotter steps (M < 6) which allows an 
accurate diagonalization of the matrices involved in the 
Trotter decomposition. At first sight, a small M limits 
the. calculation to high temperatures. However, Delica et 
al.p extrapolate their results for M = 4,5 and 6 to higher 
values of M and thus obtain thermodynamic quantities 
that are expected to be accurate to within a few percent 
in the temperature region T > 0.16 J ~ 4 K. This re- 
striction is not crucial for application to CSN1F3 in view 
of the 3D-ordering transition below Tjv = 2.7 K which 
limits the validity of the ID model anyway. 

Our calculation is based on the recently devel- 
oped transfer matrix renormalization group (TMRG) 
algorithmlrEJ which concentrates on the largest eigen- 
value of the QTM and thus addresses directly the infinite- 
chain limit. Furthermore, the number of Trotter steps 
can be chosen to be large (M ~ 160) if the resulting 
huge matrices are diagonalized by a judicious truncation 
to a finite number of important states chosen in a manner 
analogous to that employed in the earlier DMRG calcula- 
tion of ground-state properties .0 The explicit numerical 
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FIG. 2. Comparison of TMRG predictions for the tem- 
perature dependence of the zero-field transverse susceptibility 
Xx with experimental data from Ref. [l(] (solid circles). The 
dashed line corresponds to the standard set of parameters of 
Eq. ( [2,E| ) and the solid line to a lower value of the gyromag- 
netic ratio (g± = 2.18). The inset illustrates the calculated 
susceptibility for two values of anisotropy, A — 9 K (solid line) 
and A = 11 K (dotted line), which lead to virtually identical 
results. 



results discussed in the remainder of this paper were sta- 
bilized to an accuracy better than line thickness, down 
to temperature as low as T = 0.02 J ~ 0.5 K which is one 
order of magnitude lower than the lowest temperature 
reached in earlier calculations. We find that the results 
of Delica et alia are reliable, within the anticipated limits 
of accuracy, whereas the more recent elaborate semiclas- 
sical calculation of Cuccoli et al.O is not very accurate 
over the temperature region of current interest. 



A. Field perpendicular to c 

We begin with the discussion of the temperature de- 
pendence of the zero-field transverse susceptibility x± 
measured sometime ago by Dupas and Renard.Ej The 
TMRG result for the standard set of parameters ( |2.5| ) is 
depicted by a dashed line in Fig. || and is seen to system- 
atically deviate from the experimental data. The agree- 
ment with experiment for this set of parameters claimed 
by Cuccoli et al.U is due to inaccuracies in their calcu- 
lation, a point that will be made more explicit in our 
subsequent discussion of the specific heat. 

Now, the transverse susceptibility x± is found to be 
largely insensitive to the specific strength of anisotropy, 
as demonstrated in the inset of Fig. g. On the other 
hand, x± depends quadratically on the gyromagnetic ra- 
tio g± and is thus very sensitive to its specific value. It 
is then important that a reasonable agreement with the 
data is achieved for the same value n± = 2.18 obtained 
by our spinwave analysis of Sec. ||, as shown by the 
solid line in the main frame of Fig. ^|. The remaining 
systematic departure from the data observed for T < 5 
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FIG. 3. Comparison of TMRG predictions for the field 
dependence of the magnetization M at selected temperatures 
with experimental data from Ref. || (solid circles). M s is the 
saturation magnetization, and the specific choice of parame- 
ters is discussed in the text. 



K could be due to a gradual onset of 3D ordering at low 
temperatures. 

The above choice of the gyromagnetic ratio is further 
challenged by comparing, in Fig. 0, the TMRG predic- 
tion for the field dependence of the magnetization, with 
experimental data taken at selected temperatures.Q The 
specific value of A chosen in Fig. [| is not important 
because the transverse magnetization is also not partic- 
ularly sensitive to the strength of anisotropy. But the 
relative low value g± = 2.18 was again important to im- 
prove agreement with the data. Yet a significant dis- 
agreement between theory and experiment is apparent 
in Fig. ||, even at relatively high temperatures. The 
lower value g± = 2.1 employed in Ref. ^ reduces but 
does not eliminate the discrepancy. An attempt to rem- 
edy this situation by incorporating a phenomenological 
interchain interaction leads to a deterioration of the cor- 
responding theoretical prediction for the zero-field trans- 
verse susceptibility.^ 

We next discuss the specific heat C = C (T, H) which 
was measured experimentally by Ramirez and Wolf.E-3 In 
fact, most of the attention was concentrated on the excess 
specific heat 



AC = C(T, H) - C(T, 0) 



(3.1) 



viewed as a function of field H at some specified temper- 
ature T. An elementary argument based on the dilute- 
magnon approximation suggests that AC is negative and 
decreases with increasing field, because the magnon dis- 
persion discussed in Sec. || increases monotonically with 
H for all values of the crystal momentum q. Neverthe- 
less, the experiment revealed that AC rises to a posi- 
tive maximum at some field H max ~ T 2 before it be- 
gins to decrease and eventually reach negative values for 
stronger fields. A possible explanation of this unexpected 
behavior could be that the dilute-magnon approximation 
breaks down in the actual temperature range of the ex- 
periment, or "nonlinear modes" are activated in addition 




4 6 
H [kG] 

FIG. 4. Comparison of TMRG predictions for the excess 
specific heat (solid lines) with experimental data from Ref. 
|l7| (solid circles) for two typical values of temperature. The 
dashed lines depict the corresponding theoretical r esu lts of 
Ref. |^ for the same set of parameters given by Eq. (2.5). 



to magnons. Whence the beginning of a long debate con- 
cerning the possible relevance o£ gine-Gordon kinks, at 
least in some approximate sense. aU 

One of the advantages of an accurate numerical algo- 
rithm such as TMRG is that potential nonlinear effects 
are automatically taken into account. Our results for the 
excess specific heat calculated for the standard choice of 
parameters given in Eq. (2.5) are depicted in Fig. ^ 
for two characteristic values, of temperature actually em- 
ployed in the experiment .ED Inspite of the overall quali- 
tative agreement, significant quantitative differences are 
apparent in Fig. [| for both values of the temperature. We 
were thus sumrised to note that the theoretical results of 
Cuccoli et alBEl for the same set of parameters, depicted 
by dashed lines in Fig. ||, are in agreement with the data 
for the specific temperature T = 7.1 K. On the othes 
hand, our results agree with those given by Delica et al.u 
for the same set of parameters, except for some minor (a 
few percent) differences anticipated by the introductory 
remarks of this Section. As mentioned already, a simi- 
lar criticism applies to the calculation of the transverse 
susceptibility by Cuccoli et al.El. We must thus conclude 
that the semiclassical nature of their method does not 
allow a completely accurate calculation in this tempera- 
ture range and the claimed agreement with experiment 
is fortuitous. 

It is now interesting to examine whether or not the 
alternative set of parameters given in Eq. (^^) may be 
used to eliminate the observed differences. In fact, our 
results quoted in Fig. |^, together with those given in Fig. 
|] of Ref. H for yet another set of parameters, suggest 
that an accurate fit of the data is not possible for any 
reasonable choice of parameters. 

Nevertheless, the main qualitative features of the ex- 
perimental data are reproduced by the theoretical cal- 
culation. Therefore, it is important to examine further 
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FIG. 5. Comparison of TMRG predictions for the excess 
specific heat, for two different sets of parameters, with exper- 
imental data from Ref. [l?] (solid circles). 



within the ID model the mechanism by which the simple 
spinwave argument given earlier in the text is reconciled 
with a positive excess specific heat. We first consider the 
quantity 



Tln(T 3/2 C) = G + G\T + G 2 T 2 



(3.2) 



where the expansion in the right-hand side presumes that 
the low-temperature thermodynamics is dominated by 
magnons with a, q = energy gap equal to G. A deta iled 
TMRG calculation of the left-hand side of Eq. (|J) fol- 
low temperatures down to T = 0.02J reveals a behav- 
ior that is indeed consistent with the right-hand side of 
the same equation. Putting it in more practical terms, 
an extrapolation to T = using a second-degree poly- 
nomial to fit the low-temperature numerical data yields 
estimates of the magnon gap G which are in agreement 
with the direct DMRG calculation given in Table |. A 
curious fact is that the present calculation gives values 
for the gap that are even closer to the 1/n results of Table 
|[ but this may be an artifact of the specific second-order 
interpolation scheme. 

The implied normal spinwave behavior of this easy- 
plane ferromagnetic chain should be contrasted with the 
low-teroaerature anomalies discovered by Johnson and 
BonnerEj in an easy-axis ferromagnetic r chain and re- 
cently confirmed by a TMRG calculation!^] The absence 
of such anomalies in the present model reinforces the need 
for explaining the excess specific heat in simple terms. 

In the remainder of this subsection we find it conve- 
nient to work exclusivel y wi th the rationalized param- 
eters a and h of Eq. (2.2) whereas the temperature 
t = T / J is measured in units of the exchange constant 
J. The corresponding absolute specific heat per lattice 
site is denoted by c = c(r, h) and the excess specific heat 
by 8c = c(t, h) — c(r, 0Y. 

The inset of Fig. o illustrates the calculated tem- 
perature dependence of the specific heat c for a typical 
anisotropy a = 0.5 and two field values; h — and 0.1. 
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FIG. 6. TMRG calculation of the excess specific heat 
8c for a typical anisotropy (a — 0.5) and for the isotropic 
spin-1 ferromagnetic chain (a = 0). The inset depicts the 
temperature dependence of the absolute specific heat c for 
two field values, h = and 0.1, and anisotropy a — 0.5. All 
quantities shown in this figure are expressed in rationalized 
units. 



It is clear that a nonzero field causes a depression of the 
specific heat at low temperatures thanks to the open- 
ing of a finite magnon gap. This is the expected nor- 
mal spinwave behavior, as predicted by the usual dilute- 
magnon approximation. What is not accounted for by 
dilute magnons is the crossing of the h — and h = 0.1 
curves at a point P that corresponds to a specific temper- 
ature r which depends on h. In particular, P is located 
near the origin for small h and moves outward with in- 
creasing h. This crossing is precisely the origin of the 
positive excess specific heat at low h, as demonstrated 
again by the a = 0.5 solid curve in the main frame of 
Fig. |]for the specific temperature r = tq = 0.2. 

Indeed, for any fixed To, the crossing point P occurs at 
some t < To for sufficiently weak fields, and thus leads 
to positive 8c at r = tq. With increasing field the point 
P moves to the right and the corresponding temperature 
r eventually overtakes To, thus leading to negative 8c at 
t = To for sufficiently strong fields. The described picture 
is valid for any choice of To, and is confirmed by all of 
our numerical experiments. Therefore, the explanation 
of a positive 8c at low fields is equivalent to ascertaining 
the robust enhancement of the absolute specific heat c 
with increasing field, inspite of its initial depression by 
the field dependent magnon gap. 

At this point one could invoke the popular sine-Gordon 
approximation to argue that the crossing mechanism de- 
scribed in the preceding paragraph is due to the acti- 
vation of kinks or other nonlinear modes in addition to 
magnons. We think that such an interpretation is du- 
bious simply because the same mechanism occurs also 
in the isotropic Heisenberg chain, as illustrated by the 
a = line in Fig. ^. In fact, the effect is strongly pro- 
nounced in the isotropic limit, even though a sine-Gordon 
approximation is clearly out of question. 



G 



Therefore, we return to the described crossing mech- 
anism And attempt to explain it by more elementary 
means.cil The absolute specific heat satisfies the obvious 
identity 



dr c(t, h) — u(oo, h) — u(0, h), 



(3.3) 



where u(t, h) is the internal energy at temperature r and 
field h. A corresponding identity for the excess specific 
heat is obtained by applying Eq. (|3.3[) twice: 



dr 5c(t, h) = [u(oo, h) - u{oo, 0)] + [u(0, 0) - u(0, h)] . 

(3.4) 



A significant simplification occurs in the limit of an 
isotropic ferromagnetic chain for which the field depen- 
dence of the energy levels is simply a linear Zeeman shift 
mh, with m = 0, ±1, ±2, . . .. Therefore the field depen- 
dence averages out of the infinite-temperature internal 
energy w(oo,/i), which is the sum of all energy levels, 
and it(oo, h) — u(oo, 0) = 0. If we further recall that 
e(h) = u(0, h) is the ground-state energy at field h, we 
obtain the elementary sum rule 



dT8c(T,h) = e(0)-e(h) 



(3.5) 



where we have also invoked the known energy of the fully 
polarized ferromagnetic ground state. 

The obvious consequence of Eq. (3.5) is that posi- 



tive values of 5c are the rule rather than the exception. 
In particular, the initial depression of the specific heat 
(5c < 0) at low temperatures, due to the opening of a 
magnon gap at finite field, is overwhelmed by positive 
values of 5c attained at higher temperatures also thanks 
to the applied field. This explains the gross features of 
the crossing mechanism described earlier in the text and 
concludes our discussion of the excess specific heat. 



B. Field parallel to c 

The case of a field parallel to the c-axis is equally in- 
teresting but the corresponding experimental work has 
not been as extensive. We begin with the discussion of 
the temperature dependence of the zero-field longitudi- 
nal susceptibility. A notable feature of X\\(T) is that it 
must approach a finite value in the limit T — > 0. A 
simple estimate of this value is obtained by a straight- 
forward classical argument. In the presence of a field 
H = (0, 0, H) the classical ground state is such that all 
spins form an angle 9 with the c-axis calculated from 
cos 6* — g\\fiBH/2A. Therefore, the T = magnetization 
is given by M = Ngu^B cos 8 and the susceptibility by 



1 



xf(T=0) = — (N 9 p%) 



(3.6) 
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FIG. 7. Comparison of TMRG predictions for the temper- 
ature dependence of the zero-field longitudinal susceptibility 
X|| with experimental data from Ref. hq (solid circles). 

where N is the total magnetic sites and g\\ is the gyro- 
magnetic ratio for a field applied along the c-axis. 

Of course, numerical estimates based on the above clas- 
sical result are not expected to be accurate, for reasons 
similar to those explained in Sec. [n]. However, a more 
accurate prediction may again be obtained through the 
1/n expansion. To leading order, the T = magnetiza- 
tion is calculated as the expected value of the azimuthal 
spin in the Hartree variational ground state given in the 
Appendix of Ref. ^. Restricting that calculation to weak 
fields one may extract the T — longitudinal suscepti- 
bility 



l/n frr, 

X|| ( T 



0) = 



(3.7) 



The main difference from Eq. (3.6) is an overall factor of 
2, which is essentially the same factor that caused the low 
estimate A — 4.5 K in the early literature,lil in addition 
to some mild dependence on the exchange constant. In 
any case, the main conclusion is that x\\ is more sensitive 
to the value of the anisotropy constant A than to the 
exchange constant J, a situa tion that is reverse to the 



one encountered in Sec. Ill A 



Therefore, the longitudinal susceptibility is an ideal 
physical quantity to yield a sensible estimate of the 
anisotropy constant A, provided that an accurate value 
for gu is also available. The latter is fixed here by 
appealing to a theoretical estimate^!] of the difference 
g± — gn ~ 5 x 10~ 2 which leads to g\\ = 2.13 if we 
adopt our earlier value for the transverse gyromagnetic 
ratio g± — 2.18. The corresponding TMRG calculation 
of X\\{T) is illustrated in Fig. ^ for various reasonable 
choices of A. The experimental dataE3 are well repro- 
duced for the set of parameters 

J = 23.6K, A = 8.25K, gy = 2.13, (3.8) 

which is closer to the set employed by Delica et al.@. In 
addition, the field dependence of the magnetization mea- 
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H [kG] 

FIG. 8. Comparison of TMRG predictions for the field 
dependence of the magnetization M at selected temperatures 
with experimental data from Ref. ^ (solid circles). The field 
is applied along the c-axis and M s is the saturation magneti- 
zation. 

sured at selected temperatures^ agrees with our TMRG 
calculation without further fit of parameters, as demon- 
strated in Fig. @. 



Incidentally, for this choice the classical result (3.£) 
yield s 0.10 emu/mol and the leading 1/n approximation 
( |3.7|) gives 0.19 emu/mol. These values should be com- 
pared with X\\(T = 0) ~ 0.175 emu/mol extracted by a 
visual extrapolation of the solid curve in Fig. [7] to T = 0. 
Including the 1 j n co rrection produced by zero-point fluc- 
tuations in Eq. (3.7) will bring its prediction to the same 
level of accuracy with the magnon gap discussed in Table 

s 

It is now interesting to take this calculation into the re- 
gion of strong fields where the ground state becomes com- 
pletely ordered along the c-axis. Such a ferromagnetic 
state is actually an exact eigenstate of the Hamiltonian 
for any strength of the field H. But the corresponding 
magnon gap 



G = g\\HBH - A 
is positive only for H > H c where 



A/g\\HB 



(3.9) 



(3.10) 



is the critical field beyond which the fully ordered state is 
the absolute ground state. The gap vanishes for all H < 
H c because the corresponding magnon is a Goldstonc 
mode associated with the axial symmetry for this field 
orientation. 

For the set of parameters ( |3.8| ) one finds that H c = 58 
kG, in reasonable agreement with the value 62.5 kG esti- 
mated from an experiment of A. Miedan which is quoted 
in Ref. [l6] but is apparently unpublished. According to 
the description of Dupas and Renard,E3 Miedan measured 
the field dependence of the magnetization at T — 4.2 K 
and extracted H c from the observed bending of the M(H) 
curve. Although we do not know the details of this ex- 
periment, we have calculated the M(H) curve at T = 4.2 
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J = 23.6 K, A = 8.25 K 



T = 2.4K 
T = 4.2K 



0.5 



H/H c 
1.5 



H/H„ 



FIG. 9. TMRG calculation of the field dependence of the 
magnetization M for a wide field range and two typical values 
of temperature. The inset displays the corresponding results 
for the field dependence of the susceptibility. The field is 
applied along the c-axis and the critical field H c is estimated 
to be 58 kG for g« = 2.13. 



K for a wide field range and the result is depicted by a 
dashed line in Fig. [| Interestingly, the bending of the 
M{H) curve is not predicted to be especially sharp at this 
temperature, as is apparent in the corresponding suscep- 
tibility displayed also by a dashed line in the inset of Fig. 
^. In other words, if the location of the maximum of the 
susceptibility were taken as an estimate of the critical 
field H c , the latter would have been severely underesti- 
mated. The situation improves slowly at lower temper- 
atures, as indicated by the solid lines in Fig. || which 
correspond to T = 2.4 K; i.e., to a temperature that is 
already below the 3D-transition temperature Jjv = 2.7 
K. 

It is clear that we cannot go farther with our theo- 
retical arguments without explicit knowledge of detailed 
experimental data on M(H) in this field region. We thus 
conclude the discussion of magnetization with a comment 
concerning an apparent contradiction between the results 
of Fig. ^ and those given earlier in Fig. |§| for lower field 
strengths. Indeed, Fig. || suggests that the magnetiza- 
tion M(H) for any given field H decreases with increasing 
temperature, as expected, while Fig. ^ indicates that a 
relative crossing occurs between any two M(H) curves. 
The resolution of this apparent paradox lies in the fact 
that the values of temperature employed in Fig. || are all 
greater than the temperature T ~ 7.5 K, at which the 
maximum of the zero-field susceptibility of Fig. |^ occurs, 
while those of Fig. ^| are smaller. 

Finally, we discuss the specific heat in a field parallel 
to the c-axis. It appears that no measurements have been 
made for this field-orientation but could prove to be feasi- 
ble in the future.Ej Our TMRG calculation of the excess 
specific heat is illustrated in Fig. [l(^ for the two values of 
temperature employed in our preceding discussion of the 
magnetization. The characteristic double peak near the 
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FIG. 10. TMRG calculation of the excess specific heat 
for a wide field range and two typical values of temperature. 
The field is applied along the c-axis and the critical field H c 
is estimated to be 58 kG for g\\ — 2.13. 



critical field H c was anticipated by earlier workEll based 
on a classical transfer matrix calculation, on the known 
exact solution for a spin-i XY chain, as well as on an ac- 
curate numerical solution for a spin-^ XXZ chain based 
on the Bethe Ansatz. The calculated double peak is also 
a clear departure from the corresponding prediction of 
the dilute-magnon approximation^!] and could eventually 
be observed in CsNiF3 . An unfortunate feature of Fig. [H] 
is that a strongly pronounced double peak is predicted to 
occur in the low-temperature region where the ID model 
is no longer applicable. 



IV. CONCLUSION 

We have presented a more or less complete calculation 
of the dynamics and the thermodynamics associated with 

The T = dynamics is ef- 



the spin-1 Hamiltonian (1.1 



ficiently described by an Xjn expansion whose full poten- 
tial has not yet been explored. For example, an accurate 
calculation of the magnon dispersion for a field parallel 
to the c-axis is also possible but has not been carried out 
mainly because there seems to have been no experimental 
effort in that direction. 

On the other hand, the thermodynamics is calculated 
by a powerful TMRG method which has opened the 
way to obtain accurate theoretical predictions for a wide 
class of quantum magnetic chains. Suffice it to say that 
our present algorithm may be trivially adjusted to han- 
dle spin-1 Haldane-gap antiferromagnets in the presence 
of anisotropy and external fields. Even in the case of 
completely integrable spin-i chains, for which the Bethe 
Ansatz applies, the-calculation of the thermodynamics 
is far from trivial Ej Nevertheless, TMRG can be ap- 
plied in a straightforward manner irrespective of com- 
plete integrabilityEJ 



The extent to which the ID Hamiltonian (1.1) may 
describe the magnetic properties of CSN1F3 has been de- 



bated on several occasions. Our calculations confirm the 
general conclusion that the ID model accounts for the 
main features of all available experimental data. But it 
is also clear that departures from ideal model behavior 
are present, especially at low temperatures approaching 
the 3D-ordcring transition temperature TV = 2.7 K. 
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